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SUMMARY 


Simple  approximate  solutions  are  derived  for  the  relationships 
between  the  rotor  thrust  and  flight-path  velocity  components  and  the 
rotor  blade  angle,  torque,  and  in-plane  forces.  These  approximate  solu¬ 
tions,  based  upon  the  assumption  of  a  triangular  distribution  of  blade 
circulation  and  a  parabolic  variation  of  blade-element  profile  drag  with 
lift,  are  sufficiently  accurate  for  preliminary  calculations  and  the 
determination  of  the  equilibrium  angle  of  attack  and  lateral  tilt  of  the 
tip-path  plane. 

A  set  of  more  exact  blade-element  equations  is  then  derived  giving 
the  relations  between  the  thrust  and  flight-path  velocity  components  and 
the  equilibrium  blade  angles,  torque,  and  in-plane  forces  and  moments. 
Neither  the  blade-element  nor  the  approximate  solutions  are  dependent 
upon  the  usual  approximations  that  the  inflow  angle  and  blade  angle  of 
the  blade  elements  are  small  angles.  Thus  the  present  equations  should 
be  useful  for  convertaplane  as  well  as  helicopter  calculations. 

It  appears  that  nonlinear  blade  twist  may  be  desirable  for  a  con¬ 
vertaplane  rotor  in  order  to  obtain  useful  propeller  efficiencies. 
Therefore,  the  blade-element  equations  have  been  arranged  so  that  any 
reasonable  distribution  of  blade  twist  may  be  used.  Also,  the  equations 
were  set  up  in  terms  of  an  arbitrary  blade-chord  distribution  since  it 
was  found  that  the  use  of  the  actual  blade-chord  distribution  and  the 
elimination  of  the  usual  assumption  that  the  blade  airfoil  extended 
Inboard  to  the  axis  of  rotation  largely  eliminated  the  necessity  for  the 
usual  reverse-flow  corrections.  Tables  of  the  necessary  factors  are 
included  for  blades  having  a  linear  taper,  linear  twist,  and  an  airfoil 
contour  from  *r  =  0.15N  to  r  =  R  and  for  blades  having  a  linear  taper 
helical  twist,  and  airfoil  contours  extending  from  r  =  0.20R  to  r  =  R 
(where  r  is  the  radius  of  the  blade  element  and  R,  the  radius  of  the 
blade  tip) . 
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The  present  analysis  is  based  upon  the  following  assumptions: 


(l)  The  blade-element  lift  coefficient  may  be  assumed  to  be  pro¬ 
portional  to  the  sine  of  the  blade-element  angle  of  attack,  and  the 
blade-element  profile-drag  coefficient  may  be  represented  by  the  first 
three  terms  of*  a-  Fourier' series  in  the  blade-element  angle  of  attack. 
This  implies  the  neglect  of  blade  stall  effects  in  the  equations  for  the 
blade  angles.  The  effect  of  tip  stall  is  taken  into  account  in  the 
equations  for  the  rotor  torque. 


(2)  The  blade  axes  may  be  assumed  to  be,  and  to  remain,  straight  ' 


lines. 


(3)  The  lateral  and  longitudinal  variations  of  the  normal  component 
of  the  induced  velocity  at  the  tip-path  plane  may  be  assumed  to  be  linear. 


(4)  The  effects  of  compressibility  on  the  tip  sections  of  the 
advancing  blade  may  be  neglected. 


(5)  All  radial  velocity  components  and  the  tangential  components ■ 
of  the  induced  velocity  may  be  neglected. 


(6)  Blade  tip  effects  may  be  neglected, 


A  comparison  of  the  results  given  by  the  present  equations  with  the 
full-scale  helicopter  test  data  of  NACA  TN  1266  shows  good  agreement  for- 
the  helicopter  flight  range  covered  in  that  report. 


INTRODUCTION 


This  project,  which  was  conducted  at  the  Georgia  Institute  of 
Technology  Engineering  Experiment  Station  under  the  sponsorship  and  with 
the  financial  assistance  of  the  National  Advisory  Committee  for  Aeronautics 
was  undertaken  in  order  to  develop  a  blade-element  analysis  for  lifting 
rotors^ that  would  be  useful  for  convertaplane  calculations.  This'  neces¬ 
sitated  the  elimination  of  the  usual  approximations  that  the  blade-element 
inflow  angle  0  and  the  blade  angle  6  are  small  angles  and  required  a 
reasonably  exact  treatment  of  the  blade  geometry* 


...  K  was  found  that  a  practical  approach  to  the  problem  of  eliminating 
the  small-angle  approximations  for  the  lift  forces  could  be  obtained  by 
writing  the  lift  coefficient  of  the  blade  element  as 


l  -  a  sin  ar  =  s  (sin  0y  cos  0y  +  cos  0y  sin  0y) 


m 
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Mi 
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and,  consequently,  the  thrust  component  of  force  dL  cos  0V  on  a  blade 
element  as 

dL  cos  0V  =  ~  pa  (u  cos  0v)|sin  £v(u  cos  0V)  +  cos  @v(u  sin  0V)J  c  dr 


Similarly,  the  tangential  component  of  the  lift  on  a  blade  element  may  be 
expressed  as 


dL  sin  0V  =  -  pa(u  s in ' 0V)  jjs in  6 v (u  cos  0y)  +  cos  0v(u  sin  0vjj  c  dr 


It  was  also  found  that  the  small-angle  approximations  could  be  largely 
eliminated  for  the  profile-drag  terms  by  expressing  the  blade-element 
profile-drag  coefficient  c^  as 


C(*o 


sl-n  0^+62  cos  ar 


The  exact  blade  geometry  has  been  retained  throughout  by  expressing 
the  blade-chord  and  blade-twist  distribution  in  the  form  of  the  following 
constants: 


nc 


cos  e^x 


n-1 


dx 


and 


ns 


dx 
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4, 

where  6 ^  is* the  blade  twist  in  the  angle  of  zero  lift  between  the 
reference  station  and  nondime ns ional  radius  x.  Values  of  these  con¬ 
stants  are  given  in  tables  1  to  3  for  blades  having  linear  taper,  linear 
twist,  and  X]_  =  0.15  and  in  tables  1,  4,  and  5  for  blades  having  linear 

taper,  helical  twist,  and  =  0.20. 

The  present  system  of  equations  has  been  set  up  with  respect  to  tip- 
path-plane  coordinates  or  coordinates  based  on  the  virtual  axis  of  rota¬ 
tion  (fig.  l)  rather  than  the  usual  coordinate  system  based  on  the  plane 
of. zero  feathering  in  order  to  obtain  shorter  expressions  for  the  in-plane 
rotor  forces  and  moments.  The  use  of  coordinates  alined  with  the  virtual 
axis  of  rotation  also  facilitates  the  treatment  of  some  accelerated  flight 
problems. 

Certain  refinements  in  the  induced- velocity  theory,  as  given  in 
reference  1,  have  been  incorporated  with  some  minor  changes  in  the 
present  equations  along  with  the  necessary  terms  for  an  arbitrary  angular 
velocity  of  roll  and  pitch  of  the  tip-path  plane. 

Standard  NACA  nomenclature  has  been  used  where  possible,  with  the 
subscript  v  for  virtual  axis  of  rotation  appended  to  the  usual  symbols 
which,  in  this  paper,  have  a  similar  meaning  but  different  numerical 
values. 


NOTATION 

a.  slope  of  lift  curve  for  blade  element  at  0.75R  (per  radian) 

a 0  rotor  coning  angle 

a0  coning  angle  for  zero  blade-root  bending  moment 

coefficient  of  sine  component  of  blade  cyclic-pitch  angle 
measured  with  respect  to  tip-path  plane  where 

6y  =  Aq  +  6t  -  a-j_  sin  \|r  +  b^  cos  ^ 

also  coefficient  of  cosine  term  of  Fourier  series  for 
blade  flapping  angle  p  measured  with  respect  to  plane 
of  zero  feathering  where 

P  =  aQ  -  a^  cos  f  sin  \|f  -  a2  cos  2\(r  -  b2  sin  2^  -  .  .  . 
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*» 


a2 


% 

Ao 


coefficient  of  second-harmonic  cosine  term  in  Fourier  series 
for  blade  flapping  angle 

mean  blade  pitch  angle  at  reference  station,  positive  above 
tip-path  plane 


b  number  of  blades  in  rotor 


<* 


b^  coefficient  of  cosine  component  of  blade  cyclic-pitch  angle 

measured  with  respect  to  the  tip-path  plane;  also  coeffi¬ 
cient  of  sine  term  of  Fourier  series  for  blade  flapping 
angle  measured  with  respect  to  plane  of  zero  feathering 

b£  coefficient  ,of  second- harmonic  sine  term  in  Fourier  series 

for  blade  flapping  angle 

c  blade  chord  at  radius  r 


cdo 

cl 


Cmy 


extended,  blade-root  chord  at  r  =  0  (for  linear  taper) 
section  profile-drag  coefficient 
section  lift  coefficient 

rotor  rolling -moment  coefficient  measured  about  X-axis 

rotor  pitching-moment  coefficient  measured  about  Y-axis 


CQ 

acQs 


Q 


rotor  torque  coefficient  . 
increment  to  Cq  from  tip  stall  on  retreating  blade 


rotor  thrust  coefficient 


T 


PrtQ2R4 


'X 


rotor  X-force  coefficient 


x 


i  p  k 
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rotor-blade  tangential-force  coefficient,  positive  in 


direction  of  rotation 


xy 


rotor  Y-force  coefficient 


rotor-blade  thrust -force  coefficient 


fuselage  and  wing  drag 
blade  profile  drag 


mean  blade  drag  angle,  positive  in  direction  of  rotation  and 
measured  between  blade  axis  and  line  through  rotor  axis 
of  rotation  and  drag  hinge  (i.e.,  blade  drag  angle  £  is 
i  =  Eq  +  Ei  cos  +  Fi  sin  if  +  .  .  . ) 

coefficient  of  cosine  term  in  expression  for  blade  drag 
angle 

coefficient  of  sine  term  in  expression  for  blade  drag  angle 

component  of  rotor  resultant  force  acting  along  X-axis 

tangential  component  of  the  resultant  air  force  on  blade, 
positive  in  direction  of  rotation 

component  of  rotor  resultant  force  acting  along  Y-axis 

Z  component  of  resultant  air  force  on  blade 

acceleration  due  to  gravity 

mass  moment  of  inertia  of  blade  about  flapping  hinge 
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H 

k 


a0 


Lp 

Mx 

My 

Q 

r 

r 

rl 

r3 

R 


mass  moment  of  inertia  of  rotor  about  virtual  axis  of 
rotation 

mass  moment  of  inertia  of  blade  about  drag  hinge 

blade-root  spring  constant  (blade-root  bending  moment  in 
foot-pounds  divided  by  angular  deflection  in  radians  of 
.  three- quarter-radius  point  from  Sq) 

fuselage  and  wing  lift 

rotor  rolling  moment  ■ 

rotor  pitching  moment 

rotor  torque,  negative  in  direction  of  rotation 

radius  of  blade  element  c  dr 

v 

radius  of  blade  center  of  gravity 
radius  of  inboard  blade  airfoil  element 
radius  of  flapping  hinge 
radius  of  blade  tip 


t  =  y.1.?.  _  1  (for  linearly  tapered  blades) 


U 


rotor  thrust,  component  of  rotor  resultant  force  along 
Z-axis 

component  of  resultant  velocity  at  blade  element  that  is 
normal  to  blade  axis 

mean  normal  component  of  induced  velocity  at  tip-path 
plane  (positive  down  and  to  rear) 


velocity  along  flight  path 

Z  component  of  induced  velocity  at  radius  r  and  azimuth 
angle  i|r  (positive  in  plus  Z-direction) 

slope  of  longitudinal  variation  of  nondimens ional  induced 
velocity 
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gross  weight  plus  down  component  of  any  acceleration  force 
acting  on  aircraft 

i 

nondimensional  blade  radius  (r/R) 

nondimens ional  radius  outboard  of  which  retreating  blade  is 
stalled 

nondimensional  radius  of  inboard  blade  airfoil  element 

slope  of  lateral  variation  of  nondimensional  induced  velocity 

angle  of  attack  of  fuselage  measured  between  flight-path 
velocity  vector  and  longitudinal  fuselage  axis 

blade-element  angle  of  attack  measured  from  line  of  zero 
lift 

angle  of  attack  of  tip-path  plane  measured  in  the  XZ-plane 
between  flight-path  velocity  vector  and  tip-path  plane, 
positive  below  tip-path  plane 

blade  flapping  angle  at  azimuth  angle  (for  tip-path  plane, 
(3V  =  a0  -  a2  cos  2i|r  -  b2  sin  2i|r  -  .  .  . ; 
for  plane  of  zero  feathering, 

=  a0  ~  al  cos  ^  ~  ^1  sin  ♦  “  a2  cos  ^  "  b2  sin  2^  -  .  .  .) 

circulation  of  blade  element  at  radius  r  and  azimuth 
angle  \|r 

constants  in  expression  for  P  where  P  =  (rQ  +  r1  sin  i|r)x 

value  of  c^Q  at  c^  =  0 

constant  in  power  equation  for  c<i0 
(i  • e  • ,  cd0  =  6q  + 

constants  for  first  three  terms  of  Fourier  series  expressing 
relation  between  c^Q  and 

(i.e.,  c^0  =  €q  +  sin  ar  +  €2  cos 

or  C£0  -  sin  ar  +  e2  cos 

blade  drag  angle  at  azimuth  angle  ty,  positive  in  direction 
of  rotation 


2N 


9 


* 


* 


w 
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St 


twist  in  zero-lift  chord  line  between  axis  of  rotation  and 
blade  tip  for  blades  with  linear  twist,  positive  for 
increased  angle  at  tip  (i.e.,  0^  =  9jx ) 

twist  in  rotor  blade  angle  of  zero  lift  between  reference 
station  and  radius  r,  positive  for  larger  angle  outboard 

design  helix  angle  at  tip  of  blade  for  blades  with  helical 

..  twist 

pitch  angle  of  blade  element  at  radius  r  and  azimuth 
angle  measured  between  zero-lift  chord  line  and  tip- 
path  plane,  positive  above  tip-path  plane 
(A0  +  -  ap  sin  t  +  b-]_  cos  ✓  ' 

angular  displacement  of  tip-path  plane  about  X-axis  from 
horizontal 


ey 

Xv 

M-v 

P 


angular  displacement  of  tipTpath  plane  about  Y-axis  from 
horizontal 

inflow  velocity  ratio  at  center  of  tip-path  plane 
sin  <Xy  -  v1 


(V  sin  CLy  -  v\ 

~  Hr  / 


/  V  cos  Ov  \ 

in-plane,  velocity  ratio  at  tip-path  plane  .  I - ] 

\  HR  J 

density  of  air 


constants  which  express  blade-chord  distribution 

/  1  P1 

i.e. ,  =  —  /  c  dx 


itR 


X1 


» 
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c  cos  dx 


> 


sin  G-fcX11-!  dx 


constants  which  express  blade -chord  and 
blade-twist  distribution 


0c 

K 


angle  between  flight  path  and  horizontal,- positive  below 
horizontal 

inflow  angle  at  blade  element  measured  in  plane  perpen¬ 
dicular  to  blade  axis  and  between  tip-path  plane  and 
relative  wind,  positive  below  tip-path  plane 

azimuth  angle  of  blade  axis  measured  about  Z-axis  from  X-axis 
(This  angle  is  very  nearly  but  not  identically  equal  to  the 
equivalent  angle  in  plane  of  zero  feathering. ) 

ratio  of  angular  velocity  of  roll  of  tip-path  plane  about 
X-axis  to,  n 


ratio  of  angular  velocity  of  pitch  of  tip-path  plane  about 
Y-axis  to  £2 


0, 


mean  angular  velocity  of  rotor  blade  axis  about  Z-axis 


All  angles  are  in  radian  measure. 


ANALYSIS 


Value  of  Normal  Component  of  Induced  Velocity  at  Radius  r 

and  Azimuth  Angle  >| r 


It  is  shown  in  reference  1  that  for  a  lightly  loaded  single  rotor 
composed  of  a  large  number  of  blades  b  each  having  a  circulation  given 
by  the  expression 


r  =  r0  +  sin  y 


(1) 


11 


NACA  TN  2656 


the  mean  value  of  the  normal  component  of  the  induced  velocity  is 


|  arct 

(1  -  |  nv2)\|xv2  nv2 


1  -  I 


in  the  denominator  of 


Equation  (2)  was  derived  on  the  assumption  that  the  wake  extended 
to  infinity  and  had  the  form  of  a  straight  elliptic  cylinder.  Thus,  for 
those  flight  conditions  where  a  "vortex  ring”  type  flow  exists,  equa¬ 
tion  (2)  is  not  applicable  and  the  va'lue  of  v  must,  at  present,  be 

obtained  from  experiment.  The  term  ~  in  the  denominator  of 

equation  (2)  arises  from  the  lateral  dissymmetry  in  the  blade  circulation 
that  is  required  for  rolling-moment  equilibrium,  and  this  term  is  the 
only  correction  which  the  elementary  theory  makes  in  Glauert’s  original 
hypothesis  that  v  =  T/2pAV!,  where  V1  is t the  resultant  velocity  at 
the  center  of  the  rotor. 


If  the  distribution  of  the  normal  component  of  the  induced  velocity 
Vi  over  the  tip-path  plane  is  denoted  by  a  power  series  in  the  non- 
dimensional  radius  x  and  a  Fourier  series  in  the  azimuth  angle  ^  such 
that  for  the  first-order  terms 


Vi  v  ,  .  . 

—  = - +  wx  cos  f  +  yx  sin  f 

mm  J 


it  can  be  shown  from  the  results  of  reference  1  that 


»-!|c 


1  -  1.8^. 


M2 

Pv) 


\y\2  y 


y  ~  2m,v  fiR 


For  level  flight  and  (j.  >  0.15  the  expression  for  y  may  be  simplified 


y  ~  Cr 
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It  may  be  noted  that  for  a  pair  of  equally  loaded,  coaxial,  counter¬ 
rotating  rotors,  the  values  of  w  and  y  are 


y  =  0 


(8) 


Approximate  Values  of  Rotor  Blade  Angles,  Torque,  and 
X-Force  and  Y-Force  Coefficients 

It  is  convenient  for  preliminary  calculations  and  checking  and  nec¬ 
essary,  in  the  general  .case,  for  the  determination  of  the  angle  of  attack 
and  lateral  tilt  of  the  tip-path  plane  to  have  simple  expressions  of  use¬ 
ful  accuracy  for  the "rotor  torque,  X  force,  and  Y  force  that  are 
independent  of  the  rotor  blade  angles.  One  such  set  of  equations  which 
take  into  account  all  the  principal  variables  including  the  primary 
effects  of  the  reverse-flow  region  may  be  obtained  from  a  consideration 
of  the  distribution  of  the  blade  circulation.  It  may  be  noted  before 
proceeding  that  the  use  of  any  radial  blade-circulation  distribution 
other  than  the  uniform  value  assumed  in  the  derivation  of  the  induced- 
velocity  equations  will  underestimate  the  induced  torque.  Thus,  it  is 
theoretically  incorrect  to  calculate  the  induced  torque  from  blade- 
element  equations.  However,  for  the  extreme  case  of  a  triangular  dis¬ 
tribution  of  circulation  along  the  radius  and  nv  =  0.5,  the  error  in 

the  induced  torque  is  only  3^  percent  and  thus  is  probably  within  the 

errors  of  present  equations  for  the  induced  velocity. 

For  the  present  purposes  a  triangular  distribution  of  blade  circu¬ 
lation  along  the  radius  and  a  sinusoidal  variation  with  azimuth  angle 
is  sufficiently  accurate  and  will  be  used.  Then 


r  *  (r0  +  r1  sin  ^)x 


(9) 
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Writing 


U  cos  0V 


=  &R(x  +  fiv  sin  \|r) 


and 

U  sin  =  £}r|\v  +  yx  sin  \|r  +  (wx  -  aQM.v)cos  i|r 
it  follows  for  thrust  and  rolling-moment  equilibrium  that 


(10) 


(11) 


!  3jt£}R2Ct  /  .  \ 

r  - - (1  -  ~  sin  ^)x  ( 12) 

b(l  -  dv2)  V  3  J 

The  -  reference  blade  angle  Aq  at  r  =  0.75R  corresponding  to  the 
average  value  of  the  circulation  and  inflow  angle  at  this  station  and 
a  weighted  chord  may  be  obtained  from  the  substitution 

2lTao .  75R  =  ^c^O-75R  =  2(r/cU)0>^R 

or 


Aq  =  sin“^- 


Ct  j  +  16 

jtbc^l  -  nv2)  ^ 


Hv 


k  +  2^60 
729 


(13) 


where 


cxn“l  dx 


m 


For  linearly  tapered  blades  the  values  of  crn  may  be  obtained  by 
interpolation  from  table  1. 
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The  values  of  and  b]_ 

circulation  and  inflow  angle  at 
for  a^  and  i|r  =  0  and  i|r  =  rt 


obtained  from  the  differences  in  blade 

r  =  0 . 75R  for  i|r  =  jt/2  and  \|r  =  3*/2 
for  b^  are 


al 


y  -  y  8cTnv 

(l  -  |  ,vj  3*ba3  (l  -  22  ^v2  +  £  nv^ 


(15) 


and 


bi  =  -w  +  -  aoHv 


(16) 


where  the  value  of  3q  for  blade  root  moment  equilibrium  on  a  blade 
with  the  flapping  hinge  at  the  axis  of  rotation  is  approximately 


bq  = 


4bi!(i  -  nv2) 


(17) 


Similarly 


^QOM-v^ 

9  -  8nv2 


(18) 


and 


b2  =  0 


(19) 


The  value  of  the  blade-element  profile-drag  coefficient  may  be 
represented  with  sufficient  accuracy  for  the  present  purposes  by  two 
terms  in  a  power  series  in  the  blade-element  lift  coefficient  cj  such 
that 


cdo  =  B0  +  €CZ2 


(20) 
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where  for  conventional  airfoils  Bq  a  0.0080  and  e  «  O.OO80.  Making 
the  necessary  substitutions  and  integrations  the  value  of  the  rotor 
torque  coefficient  is 

cT\Jr  y Mv  -  xJ  ,  r  . 

Cq  =  — — — — — —  +  —  b&o  a4  +  —  +  [i-v^)a2 

1  -  Hv2  L 


+  ACqs  from  equation  (89b)  (where  applicable) 

(21) 


Similarly,  neglecting  blade-shank  drag  which  is  assumed  to  be  included 
in  the  helicopter  parasite  drag,  the  values  of  the  X-  and  Y-force 
coefficients  are 


/2C\2  (l  *  |  Uvg)t»4 
2W  (1  -  Mv2)2 
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angle  are  sufficiently  accurate  for  preliminary  calculations,  checking, 
and.  the  determination  of  the  angle  of  attack  and  lateral  tilt  of  the 
tip-path  plane  provided  there  are  no  large  areas  of  the  rotor  outside 
the  reverse -flow  region  that  have  blade  elements  operating  in  the  stalled 
or  negative  thrust  range.  This  implies  a  reasonable  blade  twist  for  the 
flight  conditions. 

Table  6  shows  a  comparison  of  the  values  of  the  parameters 
calculated  from  the  above  circulation  equations  with  the  flight  test 
results  of  reference  2. 


Determination  of  Angle  of  Attack  and  Lateral  Tilt 
of  Tip-Path  Plane 

Given  the  values  of  the  flight-path  velocity  V,  climb  angle  0C, 

gross  weight  and  vertical  component  of  the  inertia  force  W,  fuselage 
and  wing  drag,  lift,  moment  characteristics,  and  position  of  center  of 
gravity,  the  fuselage  angle  of  attack  and  thus  the  fuselage  and  wing 
lift  Lf  and  drag  Dp  can  be  obtained  for  the  trim  condition  by 
setting  the  summation  of  moments,  acting  on  the  fuselage  and  wing  and 
taken  about  the  rotor  hub,  equal  to  zero.  Since  the  lateral  tilt  of 
the  tip-path  plane  has  a  negligible  effect,  it  follows  from  the  geometry 
of  the  above  forces,  as  shown  in  figure  2,  that 


Dp  cos  0C  -  Lp  sin  0C  +  Fx  cos  0 

tan  6y  = - — —  (24) 

W  -  Lp  cos  0c  “  I>F  sin  0c  +  Fx  sin  0y 


is  a  good  approximation  for  unaccelerated  flight.  It  may  be  noted  that 
Dp  should  include  an  allowance  for  rotor-hub  and  blade -shank  drag.  In 

general,  the  terms  involving  Fx  will  have  only  a  small  effect  on  the 
value  of  0y  and  a  sufficiently  exact  solution  can  be  obtained  on  the 
second  iteration.  Thus,  as  a  first  approximation, 
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tan  0y  =  - 

W  -  Lp  cos  0C  -  Dj. 


sin  <pQ 


(25) 


<*v  =  0c  +  a 


(26) 


Ciji  = 


W  -  Lp  cos  f>c  -  Dp  sin  0( 
__ 

pit2  R  cos  6y 


(27) 


= 


V  COS  Qry 

___ 


(28) 


I 


= 


V  sin  Oy 

Hr" 


v 

or 


(29) 


The  values  of  v/2R  may  he  obtained  from  equation  (2)  or  by  double 
interpolation  from  table  7  which  includes  the  experimental  values  for 
vertical  descent  from  reference  3  and  estimates  of  the  values  for  the 
inclined  flight  "vortex  ring"  states.  The  values  of  w,  y,  and  Fx 

can  then  be  determined  from  equations  (4),  (5),  and  (22),  and  from  these 
second  approximations  to  the  values  of  6y,  dy,  and  pv  can  be  made 

from  equations  (24),  (2 6),  and  (28).  If  necessary,  a  new  value  of  CT 
may  then  be  obtained  from  the  equation 


* 


(30) 


and  thus  the  more  exact  value  of  Xv  from  equation  (29). 

For  helicopter  calculations  the  first  approximation  for  Ctp  is 
sufficiently  accurate,  and  if  |iv  is  small  (i.e.,  <  0.15)  the  effect 

of  Fx  on  av  may  be  neglected  for  level  flight. 
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The  tail-rotor  thrust  TT 
main  rotor  is 


required  for  a  helicopter  with  a  single 


Tt  =  Q/z 


(31) 


where  Z  is  the  perpendicular  distance  between  axis  of  main  and  tail 
rotors  and  the  value  of  Cq  may  be  obtained  from  equation  (21) .  The 

lateral  tilt  0X  of  the  tip-path  plane  for  a  single-rotor  aircraft  in 
unaccelerated  flight  is  thus 


1  C 

2  y 


CQ  T 


X  ~ 


Cip 


where  Cy  is  given  by  equation  (23). 

Application  of  Two-Dimensional  Airfoil  Theory  and  Data 
to  Rotor- Blade -Element  Calculations 
Two-dimensional  thin-airfoil  theory  demonstrates  that 


(32) 


n 


9 


*  ' 


CZ  =  a  sin  a  (33) 

For  a  two-dimensional  cascade  of  airfoils,  equation  (33)  is  modified 
by  a  multiplying  function  of  the  solidity,  chord  spacing,  and  blade  angles 
that  is  very  nearly  unity  for  average  lifting-rotor  configurations  as 
shown  in  reference  k.  Thus,  within  the  approximation  that  the  radial 
components  of  flow  may  be  neglected,  equation  (33)  should  be  applicable 
for  blade-element  rotor  theory  over  the  unstalled  range  of  blade-element 
angles  of  attack.  Beyond  the  stall,  equation  (33)  is  somewhat  less  in 
error  than  the  usual  relation  cz  =  a.a  as  can  be  seen  from  figure  3 

which  Ms  a  plot  of  the  above  expressions  and  the  experimental  values  of 
cz  against  a  for  an  NACA  0015  airfoil.  The  use  of  equation  (33), 

rather  than  the  usual  approximation  that  c z  =  aa,  allows  the  thrust  and 

tangential  components  of  lift  on  a  blade  element  to  be  exactly  expressed, 
within  the  approximations  involved  in  neglecting  radial  components  of 
the  flow,  in  terms  of  the  easily- integrated  in-plane  and  normal  components 
of  the  velocity  at  the  blade  element  U  cos  0V  and  U  sin  0V.  Thus  the 

usual  approximation  that  the  inflow  angle  0V  is  a  small  angle  may  be 
eliminated.  This  may  be  demonstrated  as  follows: 


* 
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Omitting  the  negligible  component  of  the  profile  drag,  the  thrust  dT 
on  a  blade  element  c  dr  is 

dT  =  ^  pU^cc^  cos  0y  dr  (3^) 

or  since 


cz  =  a  sin  0^.  =  a (sin  0y  cos  0y  +  cos  0y  sin  0y)  (35) 


dT  ='  —  pac  (u  cos  K)  [sln  0y(u  cos  0V)  +  cos  ey(u  sin  0V)J  dr  (36) 

The  tangential  component  of  the  lift  on  a,  blade  element  may  be  similarly 
expressed  as 


dL  sin  0V  =  pac(u  sin  Is  in  0v(u  cos  0V)  +  cos  @v(u  sin  0yH  dr  (37) 


The  value  of  the  slope  of  the  lift  curve  a  of  the  blade-element  air¬ 
foil  in  the  above  relations  may  be  taken  as  the  value  corresponding  to 
the  Reynolds  number,  Mach  number,  and  surface  roughness  existing  at  the 
three-quarter-radius  point  of  the  rotor  blades  under  consideration.  For 
the  usual  tip  speeds,  in  the  500-foot-per- second  range,  the  Prandtl- 
Glauert  Mach  number  correction 


where 

af 


a ! 

^1  -  M2 


(38) 


low  Mach  number  lift-curve  slope  from  two-dimensional 
wind-tunnel  tests 


M 


free-stream  Mach  number  at  three-quarter  blade 
radius 


may  be  used  to  correct  the  lift-curve  slope  from  low  Mach  number  data. 

The  values  of  c^  obtained  from  two-dimensional  wind-tunnel  tests 
at  appropriate  Reynolds  numbers  and  model  surface  roughness  should  be 
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directly  applicable  to  rotor-blade-element  calculations  in  the  unstalled 
range  of  angles  of  attack  below  the  Mach  numbers  and  angles  of  attack 
for  drag  divergence,  since  the  effect  of  subsonic  Mach  number  on  profile 
drag  is  negligible  as  shown  in  reference  5*  However,  it  should  be  noted 
that  the  profile-drag  coefficient  is  only  constant  with  change  of  sub¬ 
sonic  Mach  number  if  it  is  taken  as  a  function  of  the  lift  coefficient. 

If  the  airfoil  section  data  are  plotted  against  section  angle  of  attack 
and  the  Prandtl-Glauert  correction  is  applied  to  the  lift-curve  slope, 
this  is  equivalent  to  multiplying  the  section-angle-of-attack  scale  by 
o  1/2 

(l  -  M2)  '  .  Consequently,  the  section-angle-of-attack  scale  on  the 

profile-drag  curve  must  be  multiplied  by  (l  -  M2)"^2  to  retain  the 
same  relation  between  c^  and  c^. 

In  view  of  the  errors  in  the  magnitude  and  distribution  of  the  blade 
circulation  that  arise  from  the  necessary  neglect  of  blade  deflections, 
and  so  forth,  it  is  probably  -not  justifiable  to  take  into  account  sec¬ 
ondary  effects  of  the  profile  drag.  Thus,  expressing  the  relation  between 
the  profile-drag  coefficient  and  the  blade-element  angle  of  attack  by  the 
first  three  terms  of  a  Fourier  series  gives 


%  =  e0  +  €1  sin  ar  +  c2  cos  ar  (39) 

The  constants  in  the  above  equation  may  be  evaluated  from  the  two- 
dimensional  wind-tunnel  data  for  the  blade  airfoil  at,  say,  a  =  0°,  5°, 
and  10°.  The  advantages  of  equation  (39)  over  the  usual  expression 


cdQ  =  s0  +  &lar  +  52ar2 

are:  The  last  two  terms  of  equation  (39)  can  be  exactly  expressed  in 
the  known  velocity  components  U  cos  0V  and  U  sin  0V;  the  resulting 

expressions  for  the  forces  and  moments  on  the  blade  are  considerably 
simplified  by  the  absence  of  the  squared  term  in  Oj.;  and  it  is  an 

equally  accurate  approximation  to  the  experimental  values  of  as 

may  be  seen  from  figure  b.  However,  in  using  equation  (39)  it  may  be 
noted  that  the  calculated  value  of  c^  is  the  small  difference  between 

large  quantities  and  thus  the  values  of  €q,  and  e2  should  be 

determined  to  four  places  in  order  to  obtain  the  value  of  c^  to  the 
customary  accuracy.  For  the  more  severe  convertaplane  flight  conditions 
where  the  inflow  velocity  is  large  (|Xv|  >  O.io)  a  certain  error  arises 
in  the  treatment  of  the  e0  terms,  and  it  is  necessary  to  fall  back  on 
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the  two-term  approximation  for  ca0,  =  e-^  sin  ctj.  +  e2  cos  ar, 

where  and  €2  are  evaluated  from  the  experimental  data  at,  say,  ■ 

ct  =  2°  and  a  =  J°.  This  additional  approximation  is  permissible  for 
these  flight  conditions,  since  the  relative  effects  of  the  profile  drag 
become  less  important  as  the  inflow  velocities  and  rotor  blade  angles 
increase.  For  example,  in  propeller  calculations  the  single-point 
approximation  c<jo  =  ec^  is  usually  used. 

It  follows  from  the  geometry  and  equations  (35)  and  (39)  that  the 
tangential  component  of  the  profile  drag  on  a  blade  element  may  be 
expressed  as 


dD0  cos 


0V  =  |  Pc(u  cos  0v){cOu  + 

el[(u  cos  0v)  s^-n  0v  +  (u  sin  0V)  cos  ®v]  + 

e2  [(u  cos  0V)  cos  6V  -  (u  sin  0V)  sin  ey]J‘  dr  (4o) 


Thrust  of  a  Blade  at  Azimuth  Angle  i|r 
The  thrust  Fz  of  a  blade  at  azimuth  angle  i|r  is 


1  fR 

Fz  =  ^  Pa  /  c 


(u  cos  0V)[(U  cos  0v)  sin  ev  +  (u  sin  0v)  cos  0v]  dr 


m 


where  is  the  radius  of  the  inboard  blade  airfoil  element.  In  the 

general  case  it  follows  from  the  geometry  that 

U  cos  0V  =12R  (x  +  4V  sin  f)  (42) 


U  sin  0V  =  fiR  (wx  +  cuyx  -  a0|iv)  cos  i|r  +  (y  -  a4c)x  sln  t  + 

2b2X  cos  21lf  -  2a2*  sin  2ij  (43) 


and 
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where 

a**  ratio  of  angular  velocity  of  roll  of  tip-path  plane  about 

X-axis  to  fl 

toy  ratio  of  angular  velocity  of  pitch  of  tip-path  plane  about 

Y-axis  to  H 

Neglecting  the  higher  harmonics  of  the  cyclic  pitch  that  may  arise  from 
control- system  linkages,  the  pitch  angle  0Y  of  a  blade  element  at 
radius  r  and  aximuth  angle  >|r,  measured  with  respect  to  the  tip-path 
plane,  is 

Gy  =  Aq  +  0^  -  a*]_  sin  cos  \|r  (44) 


where 

Aq  mean  blade  pitch  angle  at  reference  station 

0 t  twist  in  rotor  blade  angle  of  zero  lift  between  reference 

station  and  radius  r 


minus  the  coefficient  of  sine  component  of  blade  cyclic- 
pitch  angle  measured  with  respect  to  tip-path  plane 

b]_  coefficient  of  cosine  component  of  cyclic-pitch  angle 

measured  with  respect  to  tip-path  plane 

In  the  general  case  (i.e.,  for  the  convertaplane)  Aq  and  6 ^ 

may  not  be  small  angles*.  However,  it  appears  that  the  magnitude  of  the 
cyclic-pitch  angle  will  always  be  limited  by  tip  stall  on  the  retreating 
blade  to  the  range  where  it  is  a  good  approximation  that 


sin  sin  +  b^  cos 


*) 


=  -a^  sin  +  b^  cos  \|r 


(45) 


cos 


^-a^  sin  ^  +  b^  cos 


(46) 


and 
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It  follows  from  equation  (44),  upon  expanding  the  functions  sin  0, 


and  cos  0y,  that 


sin  0V  =  [sin  A0  +  (cos  A^-a^  sin  \(r  +  bx  cos  ty) j  cos  0t  + 

[cos  Aq  -  (sin  A^J^-aj  sin  \|/  +  b^^  cos  i|rjj  sin  0 ^  (47) 


cos  0V  =  [os  Aq  -  (£in  Ag)[a-^  sin  \|/  +  b-j_  cos  cos  0^  - 


[sin  A0  +  (cos  A^^aj  sin  i|f  +  b^  cos  i|/)j  sin  0 1 


(48) 


Substituting  the  values  of  U  cos  0V,  U  sin  0y,  sin  0y,  and  cos  0y 
from  equations  (42),  (43),  (47),  and  (48)  in  equation  (4l),  defining 


CTnc  =  rjr  /  cxn-1cos  6t  dx 

Jx-, 


(49) 


°ns  - 


jtR 


cxn_1sin  e^'dx 


X1 


(50) 


"^nc  CTnc  s^n  ^0  +  CTns  cos 


(51) 


■^■ns  ~  ans  sln  Aq  -  anc  cos  Aq 

multiplying  out  the  terms,  and  reducing  the  functions  of  \|r  to 

F 

harmonic  form  give  for  the  thrust  coefficient  C,  =  _ - _  of 

blade  at  an  azimuth  angle  i|f  the  expression  of  equation  (53): 


(52) 


one 


25 


NACA  TN  2656 


Equation  (53)  is  written  in  tabular  form  where  the  coefficients  in  the 
boxes  must  be  multiplied  by  row  and  column  heads.  Values  of  anc  and 
crns  may  be  obtained  by  interpolation  from  tables  2  and  3  for  linearly 
tapered  and  twisted  blades,  where  • 


c  =  cD(l  +  tx)  from  x  =  0.15  to  x  =  1  (5^-) 


et  =  elx 


(55) 


(56) 


and 

c0  extended  blade-root  chord  at  r  =  0  _ 

L  Tip  chord  -  ' 

t  =  — - -  -  1  1 

co 

6i  twist  in  angle  of  blade  zero  lift  between  axis  of  rotation 

and  tip 

In  order  to  use  the  tabulated  values  of  cxnc  and  0ns  for  blades  with 

linear  twist  and  taper,  it  is  necessary  to  take  the  reference  blade 
pitch  angle  A q  at  the  extended  blade-root  chord  cQ  at  r  =  x  =  0. 

The  use  of  the  lower  limit  x^_  =  0.15  in  the  computations  for  the 

blades  having  linear  taper  and  twist  corresponds  to  present  practice 
and  largely  eliminates  the  necessity  of  making  any  reverse-flow  correc-  t 
tion  to  the  blade  thrust.  The  reverse-flow  effects  are  discussed  in  the 
following  section. 

Additional  tables  (tables  h  and  5)  give  the  values  of  anc  and  ans 
for  blades  having  linear  taper  from  x-j_  =  0.20  to  x  =  1  and  helical 
twist  where 


et 


tan’ 


/tan 


(57) 
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and  0T  is  the  design  helix  angle  at  x  =  1.  In  this  case,  the  refer¬ 
ence  station  for  AQ  is  taken  at  the  blade  tip.  The  tables  for  helical 
twist  are  included  for  convertaplane  usage  since  helical  twist  would 
appear  to  be  desirable  for  a  reasonable  propeller  efficiency.  An  inner 
limit  of  x^  =  0.20  was  used  for  the  computation  of  the  values  of 
anc  an<^  ans  ^or  "this  case  of  helical  twist  in  order  to  minimize  the 
severe  root  stall  likely  to  occur  under  some  convertaplane  flight  con¬ 
ditions.  It  might  be  pointed  out  that  helical  twist  would  also  appear 
to  afford  an  increase  in  helicopter-rotor  performance  over  that  obtain¬ 
able  with  linear  twist. 


Reverse-Flow  Considerations 

For  normal  helicopter  and  convertaplane  flight  conditions  where 
^here  is  s  downflow  through  the  rotor  and  0y  is  negative  over  the 
reverse-flow  region,  the  maximum  value  of  ny  is  limited  for  conven¬ 
tional  rotors  to  relatively 'low  values  of  the  order  of  0.30  by  tip  stall 
on  the  retreating  blades.  Under  these  conditions  the  portion  of  the 
retreating  blade  extending  inboard  from  the  outer  edge  of  the  reverse- 
flow  region  at  x  =  sin  i|r,  where  the  in-plane  component  of  velocity 

is  zero,  to  x  =  x^,  where  the  blade  airfoil  section  ends,  has  a  neg¬ 
ligible  thrust  loading  because  the  in-plane  components  of  velocity  are 
very  small.  The  present  equations  take  into  account  the. fact  that  the 
blade  airfoil  does  not  exist  inboard  of  x  =  xj_,  for  which  region  the 
the  in-plane  components  of  velocity  are  larger,  within  the  reverse-flow 
circle,  and  previous  equations  erred  in  assuming  the  blade  airfoil  to 
exist. 

For  those  flight  conditions  where  there  is  an  upflow  through  the 
rotor  and  the  tip-stall  limitations  on  ny  are  relaxed,  the  present 

equations  give  the  proper  direction  to  the  blade-element  thrust  for 
those  blade  elements  within  the  reverse-flow  region  and  inside  the  radius 
where  0y  %  2 6y. 


Thus,  for  all  practical  purposes,  it  is  not  necessary  to  use  reverse- 
flow  corrections  when  applying  the  present  equations  to  conventional  rotors. 


For  unconventional  rotors  operating  with  net  downflow  at  large  values 
of  it  would  appear  from  strip  analysis  to  be  desirable  and  even  nec¬ 


essary  to  minimize  the  forces  in  the  reverse-flow  region  by  using  a  suf¬ 


ficiently  large  design  value  of  xj,  for  example, 


X1  >  f  ^vmax. 


In 


this  case  less  error  is  introduced  by  taking  into  account  the  inboard 

blade  airfoil  limit  and  omitting  the  usual  reverse-flow  correction  than 
vice  versa. 
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9  Mean  Rotor  Thrust 

Omitting  the  coefficients  of  the  second- harmonic  flapping  angle 
which  have  a  negligible  effect  on  the  mean  rotor  thrust,  the  value  of 
the  mean  rotor  thrust  coefficient  obtained  from  the  first  row  of  equa¬ 
tion  (53)  "by  averaging  the  value  of  Cz/a  over  the  interval  from 
ijr  =  0  to  \|f  =  2ir  is  given  by 

Ir  =  L1  +  I  a!(y  '  °*)  ■  i  bl(V  +  ^>30  +  I  a0bl^vI2c  + 


2(alNr  +  ^vVv-^lc  “  \r  al^v  +  2  (y  ~  CDx)^'v  *^2s  (58) 


Mean  Rotor  Air  Rolling  Moment 

The  value  of  the  mean  rotor  air-rolling-moment  coefficient  about 
the  X-axis 


mx  1  O  c: 
i  PrtQrR? 

is  founds  upon  integration,  to  he  obtained  by  multiplying  the  second  row 
of  equation  (53)  "by  ^  b  and  introducing  the  moment  arm  by  writing  the 
subscripts  of  Inc  and  Ins  to  one  higher  order.  Thus, 

\ 

=  |f^v  +  alxv  +  f  al(y  '  “xK  -  l  bl(w  +  *vK>3c  + 

—  aQb-j_pv  l2C  +  (sq  y  +  u>x)l4s  +  8mv  _  \t)pv^2s  (59) 

Mean  Rotor  Air  Pitching  Moment 
Similarly,  the  mean  rotor  air -pitching-moment  coefficient 

My 

^  i  p*!2i<5 
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obtained  from  the  third  row  of  equation  (53)  is 

=  [j>l*v  "  |  al(w  +  “jr)t*v  +  |  bl(y  -  %)^v]I3c  + 
l  a0ai^v21 2c  +  (bl  +  v  +  “Y^ls  "  a0^vI3s  +  \  Vv2;c2s  '  (6°) 


Mean  Blade-Root  Air  Moment 

The  coefficient  Cm0  of  the  blade-root  air  moment  Mq  is  merely 

the  first  row  of  equation  (53)  with  the  I  factors  to  one  higher  sub¬ 
script.  Thus,  for 


Mo 


'mo 


|  p*fl2R5 


Cmo 

a 


=  [i  +  |  ai(y  -  |  bi(w  +  -yj]  lie  +  I  aobi^vI3c  + 


2  (al^v  +  ^v)^v^2c  "  [>  “  al^v  +  2^  CBx)'1vJ 


3s 


(61) 


l 

Equilibrium  Values  of  Mean  Rotor  Pitching  Moment 
and  Rolling  Moment 

If  an  external  moment  M^_  is  applied  to  a  single  rotor  with  three 
or  more  blades  about  a  diameter,  axis  1,  the  differential  equations  of 
motion  about  axis  1  at  t  =  and  axis  2  at  ^  =  +i  +  90°  can  be  shown 
by  the  use  of  Euler's  equations  to  be 


1  ^ 

2  dt 


+  fiajg  + 


Mj 

Iv  ly 


(62) 
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and 


t 


1 

2  dt 


20>p 


k2^2 


■v 


(63) 


where  o)p  anc^  are  the  angular  velocities  of  the  tip-path  plane 

about  axes  1  and  2 ,  respectively,  kpGOp  and  are  the  damping 

moments,  and  Iy  is  the  mass  moment  of  inertia  of  the' rotor  about  the 

virtual  axis  of  rotation.  The  general  solution  of  equations  (62) 
and  (63)  is.  a  pair  of  equations  of  the  form 


In  the  actual  case,  damping  of  the  nutation  appears  to  be  very  rapid 
for  an  articulated  rotor.  Also,  for  pilot-controlled  motion,  ^  0. 

For  example,  for  a  constant  control  moment  Mp,  k2  =  0,  and  kp  =  2CIV, 
which  is  then  the  value  of  kp  for  critical  damping,  it  follows 'that 


(65) 

*v 

or 

M1  L  -2fit\ 

o>2  =  — -ll  -  e  )  (66) 

J-yi  l 

It  can  be  seen  from  equations  (65)  and  (66)  that  the  transients  decay 
very  rapidly  and  their  effects  can  be  neglected  in  most  problems.  There¬ 
fore,  to  a  good  approximation  for  a  single  rotor 


2 

n 


“xf 


(67) 


My  —  —  IyH  +  Myp* 


(68) 
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where  Mxf  and  My^  are  any  moments  transmitted  about  the  X-  and  Y-axes 
from  the  fuselage  to  the  rotor.  For  steady  straight  and  level  flight 


=  0  (69) 

For  steady  banked  turns  the  value  of  CT  can  be  taken  proportional  to 
sec  0X.  Also 


-g  sin  6 y  tan  Qx 

%  ~  -  (TO) 

Vft 

and 

g  sin  9X  tan  6V 

“V  ~ - m -  <71) 


where  Gx  is  the  equilibrium  lateral-tilt  angle  of  the  tip-path  plane 
(approximately  equal  to  equilibrium  angle  of  bank,  positive  for  turns 
in  direction  of  rotor  rotation). 


For  any  curvature  of  the  flight  path,  the  nondimens ional  compon¬ 
ents  and  ay  of  the  spatial  angular  velocity  of  the  aircraft  may 

be  calculated  and,  consequently,  the  approximate  equilibrium  values  of 
Mx  and  My  can  be  obtained  from  equations  (67)  and  (68). 


* 


? 


w 


Approximate  Solution  for  Equilibrium  Values  of  Mean  Reference 
Blade  Angle  Aq,  Lateral  and  Longitudinal  Components  of 
Cyclic  Pitch  a^  and  b]_,  and  Coning  Angle  aQ 

An  approximate  solution  of  the  set  of  four  nonlinear,  transcendental 
equations  (58) ,  (59) ,  (60),  and  (6l)  for  the  four  unknowns  Aq,  a. sq, 
and  bi  that  is  sufficiently  accurate  for  most  steady-flight  helicopter 
work  and  useful  as  a  first  trial  for  steady- flight  convertaplane  calcula¬ 
tions  may  be  obtained  as  follows:  Setting  the  small  terms  and  a*,  cay, 
and  Cnir  equal  to  zero  and  cos  Aq  =  1  in  equations  (58)  and  (59)  and 
eliminating  a^  gives 


sin  Aq 


(S1  '  ff3s  -  V2cJ 

\akc  + 

J  ^v2ff2c  -  V3s) 

1  +  Hv2(I2cl 

(2a3s  +  V2c) 

( 

f  1  2 

vff3c  +  f  alc  - 

V^s) 

(ff4c  +  l  ^v2<j2c  - 

V3s)  - 

2^v2cr2ca3c 

* 


(72) 
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Then,  from  equation  (59)  for  1%  =  ay  =  =  0 

_  2^vI3e  "  ,^4s  “  \r^v^2s  ,  . 

ai  _  ~  ~r  3  \  3  o  ~  (73) 

O  +  £  yuvj^c  +  I4s  +  l  ^v  I2s 

Let  he  the  design  coning  angle  for  the  general  case  of  semirigid 

blades  (i.e.,  coning  angle  for  zero  blade-root  bending  moment).  Let 
ka0  he  the  spring  constant  of  the  blade  for  angular  deflections  of  the 

three-quarter-radius  point  from  ey.  Then  setting  the  summation  of 
moments  about  the  blade  root  equal  to  zero  and  solving  for  aQ,  the 
coning  angle  nt  the  three-quarter-radius  point. 


where 


Mg  mass  of  blade 


r  radius  of  blade  center  of  gravity 


Ip  mass  moment  of  inertia  of  blade  about  flapping  hinge 

(or  root) 


|lf  the  blades  have  a  flapping  hinge  at  the  axis  of  rotation  aQ  =  kaQ  = 

If  the  flapping  hinge  is  located  at  radius  rp  from  the  axis  of  rotation 

rprMgfl2  \ 

a0  -  0  and  ka o  ~  - ; - r*  ■  Then,  knowing  an,  it  follows  from 

°  1  -  (rp/0.75R)  J  0 

equation  (60)  that  for 


%  -  ay  -  Cmy  -  0 

a0^vI3s  ~  wI4s 
V^c  +  I4s  +  ^  ^v2l2s 


(75) 
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For  those  steady,  unaccelerated  flight  conditions  where  cos  A q  2  1,  the 
above  solutions  are  sufficiently  accurate  and  may  he  used  to  calculate 
the  blade  loadings  and  rotor  torque,  X  force,  and  Y  force. 


"Exact"  Solution  for  A0,  ap,  and  bp  for  Accelerated  Flight 

Conditions  and  Those  Flight  Conditions  where  cos  Aq  ^  1 

A  reasonably  rapid  and  sufficiently  accurate  solution  of  the  "exact" 
equilibrium  equations  given  by  the  first  three  rows  of  equation  (53)  can 
be  obtained  by  using  an  approximate  value  for  the  coning  angle  Sq  such 

as  that  given  by  equation  ( IT)  or  (7*0  • 


Then  for  the  approximate  value  of  Aq  given  by  equation  (72)  and, 
for  example,  two  other  values  several  degrees  successively  smaller,  the 
"exact"  corresponding  values  of  ap  and  bp  can  be  determined  by 

rewriting  the  equilibrium  equations  for  the  rotor  pitching  and  rolling 
moments  in  the  form 


Aa]_  +  Bb1  =  P 


fmy 

ab 


where 

A  =  ^(w  +  o>y)^vl3C  -  |  a0nv2I2c  (77) 

B  =  -  jj;(y  -  ^x^v^c  "  4  ^v  I2s  ("^0 

c  =  M3c  +  f(y  -  ‘MMsc  +  J4s  +  f  ^v2i2s  (^9) 

D  =  _  I(w  +  (ny)dvI3c  +  |  ao^v2l2c  (8o) 

P  =  (w  +  0^)1^  -  aonvl3S  (81) 

R  =  -2hvI3c  +  (y  -  0^)143  +  VvI2s  (82^ 
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Then 


Having  computed  the  values  of  a^  and  b-^  for  each  of  the  assumed 
values  of  Aq,  the  corresponding  values  of  Cij.  may  be  found  from  the 
equation  for  the  thrust  equilibrium  where 

IV  =  I3c  +  2  Ilc  "  XvI2s  "  |(y  "  “x)^vI2s.+ 
j|(y  -  o^)l3c  +  |  VvIlc  +  M2s]s1  + 

|  aoM2c  -  |(w  +  aV)I3c|bl  (Q5) 

Then  plotting  the  values  of  2Cp>/ab,  a]_j  and  "bp  against  the  trial  value 
of  AQ,  the  "exact n  value  of  A q,  and  thus  ap  and  bp,  may  be  obtained 
from  the  plot  at  the  design  or  desired  value  of  Cp. 
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In-Plane  Component  of  Force  FXy  on  a  Blade 
at  Azimuth  Angle  t 

The  in-plane  component  of  force  in  the  direction  of  rotation 
on  a  blade  at  azimuth  angle  is  from  equations  (37)  anA  (4o) 


cos  e. 


(u  sin  0vJj  +  e2|cos  ©v(u  cos  0V) 


where 


cdo  =  e0  +  el  sin  ar  +  e2  cos  °!r 

Then,  by  (l)  substituting  the  previously  evaluated  expressions  for 
U  cos  0V,  U  sin  0V,  sin  6y,  and  cos  6V  given  by  equations  (42),  (43), 
(47),  and  (48);  (2)  neglecting  the  effects  of  second-harmonic  flapping; 
and  (3)  writing  (eoU)(U  cos  0V)  as  6q(U  cos  0v)^(u/U  cos  0V)  and 
expanding 


by  the  binomial  theorem  and  dropping  third  and  higher  terms,  the  expres¬ 
sion  for  the  constant  and  first-harmonic  terms  becomes 


xy  -  r 

2 


JOL 


(ACxy)a  "  (ACxy)?n  -  (ACXy)£i  -  (ACxy) 


(87) 
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Kotor  Torque 


The  effects  of  tip  stall  at  the  higher  values  of  pv  and  Crp/cr^ 
on  Cq  are  large  and  may  he  approximately  evaluated  for  high-speed 
flight  where  \v  is  negative,  as  follows:  The  retreating  blade  will 
he  stalled  outboard  of  the  nondimens ional  radius  (for  \v  negative) 


xs  = 


K  +  M-v  tan  l-  -  Aq 


al  -  A0t, 


(89a) 


-  A0  -  al  -  ^t; 


where  Aet  is  the  aerodynamic  blade  twist  between  the  reference  station 
and  the  tip.  Assuming  a  jump  of  0.08  in  the  value  of  cdQ  at  the  stall 

and  that  the  rotor  area  within  which  blade  stall  exists  is  a  segment  of 
minimum  radius  xs  and  symmetric  about  \|/  =  3^/2,  the  increment  ACQg 
to  Cq  due  to  tip  stall  is  approximately 


ba4 


acQs  *  -  ^v)2(l  -  x8)  n  -  xs2 

OJt 


(89b) 


(If  xs  <  l^v  or  xs  >  1  equation  (89b)  is  not  applicable  and 
^Cqs  =0.)  Then 

2CQ  /  "* 

-  -(Constant  terras  of  Cxy  with  subscripts  n  on  on,  Inc,  and 

2AC, 


Ins  increased  to  n  +  l)  + 


'Qs 


(90a) 


For  steady-state  calculations  equation  (90a)  may  be  reduced  to 


NACA  TN  2656 


2Cq  _  2Ct(|  y^v  -  V)  +  2ACQS  + 
b  b(l  -  tiv2)  b 


e°  K1  +  i  V^°k  ’  2  +  |(^v2  +  ^v2)  02]  + 

eijlltc  +  |  +  +  aiV)^vI2c  +  (al^v 

€  2  "  0I3c  "  (l  '  |  *iy)l4s  '  \  aobiJivI3s 

^  (al^v  +  ^v)  ^v-^-2  J 


■)X3b] 


(90b) 


Rotor  X  Force 

The  value  of  the  rotor  X-force  coefficient  Cx  is 

=  -^Sine  terms  of  Cxy^  (91) 

However,  the  greater  part  of  Cx  arising  from  the  lateral  variation  in 
blade  circulation  is  a  small  difference  between  large  quantities  which 
are  principally  functions  of  a]_  and  Aq.  It  follows  that  this  part  of 

Cx  is  more  accurately  obtained  from  the  circulation  equations  than  from 
the  blade-element  equations.  Thus  for  steady-state  solutions 

PC  2Ct(2X.v+iv  -  y) 

—  ft*  - - - p- —  +'  e0  ( 2^v  +  y^v)°r2  + 

b  b(  1  - 

el^ai^v  +  2^v)i2c  +  (al  “  y)*3s  +  (^  al^v  ■  + 

e2  [(al  "  y)l3c  +  alMv  “  ^v^M-v^lc  “  (^M-v  +  al^v)^2s]  (92) 

Rotor  Y  Force 

The  value  of  the  rotor  Y-force  coefficient  is 


2c 

=  Cosine  terms  of  Cxy 


(93) 
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As  in  the  expression  for  Cx  the  above  value  of  Cy  given  by  the  blade- 

element  equation  is  a'  small  difference  between  large  quantities  and  the 
result  for  steady-state  flight  is  more  accurately  obtained  from  the 
circulation  expression 


ct(w  -  I  ao^v) 


Second-Harmonic  Flapping  * 

Again  letting  kaQ  be  the  spring  constant  relating  the  blade-root 

bending  moment  in  foot-pounds  to  the  angular  deflection  in  radians  of  the 
three -quarter-radius  point  of  the  blades  from  the  unstressed  position, 
it  follows  that  the  magnitude  of  the  cosine  component  of  the  second 
harmonic  of  the  blade  flapping  angle  is 

a0  ~  J  +  KL  (qk\ 

&2  ~  1  -  KM  W) 

Similarly  the  magnitude  of  the  sine  component  is 

b2  x  L±M  (9 6) 

^  1  -  KM 


where 


^terms  not  involving  in  'the  cos  row  of  thrust  equa¬ 

tion  (53)  with  the  I  factors  changed  to  one  higher  sub- 

\  (  i-  \ 

script]  x  \— — - - 

\3I  l^2  -  kao/ 

^coefficients  of  bp  in  the  cos  2ty  row  of  thrust  equa¬ 
tion  (53)  with  the  I  factors  changed  to  one  higher 


subscript^ 


31  l^2  -  kg. 


(terms  not  involving  ap  in  the  sin  2 ty  row  of  thrust 
equation  (53)  with  the  I  factors  changed  to  one  higher 


subscr: 


^  pjtn^a 

31  ifl2  -  kao 


ko 
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M 


^coefficients  of  a2  in  the  sin  2^  row  of  thrust 
equation  (53)  with  the  I  factors  changed  to  one  higher 


subscript 


ipt) x  ^ 


For  steady-state  flight  conditions  where  =  ay  =  0  the  expres¬ 
sions  for  the  factors  J,  K,  L,  and  M  may  be  simplified  to 


i  pnan2R^  r- 

JZ^—2 -  - 

31x0?  -  kaoL 


2  (aPv  +  ^v)^v^2c  “  (al  "  2  ^vP3s  (91) 


K 


^  p«afi2R^ 


3ixn2  -  k 


a0 


(-2i0 


(98) 


rj  p«afi2R5 
3I]fi2  -  ka 


£  2  blA'vuvI2c  (bl  +  2  w)m3s  +  \  a0^v2l2s 


(99) 


M 


i  pnafi2R5 


3Ii«  -  *a0 


K>) 


(100) 


and  I-j_  is  the  mass  moment  of  inertia  of  the  blade  about  the  flapping 
hinge . 

It  may  be  noted  that  ka^  =  0  for  blades  having  a  flapping  hinge 

at  the  axis,  of  rotation.  If  the  flapping  hinge  is  located  at  radius  r^ 
then 


Ka0  ~ 


_  2 


0.T5R 
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Amplitude  of  Constant  and  First-Harmonic  Components  of 

Lag  Angles  in  Unaccelerated  Flight 

For  an  articulated  rotor  having  lag  hinges  normal  to  the  plane  of 
rotation  and  located  at  a  small  radius  e  the  equilibrium  blade  lag 
angle  Eq  is 


from  equation 


(101) 


where 


M 


is  the  mass  moment  of  the  blade  about  the  lag  hinge. 


Similarly  the  coefficients  of  the  cosine  and  sine  components  of  the 
lag  angle  are 


and 


E- 


I  P*R5E£  -  2ao*lsl£ 


M^e 


t 


\  +  2a0alsI^ 


M^e  -  1^ 


(102) 


(103) 


where  a-^  and  blg  are  the  cos  t  and  sin  ^  components  of  the  angle 

between  the  tip-path  plane  and  the  hub  plane..  For  unaccelerated  flight 
the  values  of  a^s  and  b-j_s  a, re  approximately 


als  *  av  *  af 


(164) 


^*ls  ~  ®xf 


(105) 


where  is  the  equilibrium  lateral  tilt  of  the  fuselage 
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Also 

h 


E£ 


mass  moment  of  inertia  of  a  blade  about  lag  hinge 

coefficient  of  cos  in  equation  (88)  for  CXy  with  sub¬ 

scripts  of  I  factors  changed  from  n  to  n  +  1;  approxi* 
mate  value  from  circulation  equations  is 


E  Y  = 


(f  w 


(106) 


coefficient  of  sin  \|r  in  equation  (88)  for  Cxy  with  sub¬ 
scripts  of  I  factors  changed  from  n  to  n  +  1;  approxi¬ 
mate  value  from  circulation  equations  is 


Fr 


(|  y  -  t  x^v) 


2CtIt  y  -  t  xv^vJ  0.0850^1^0^ 

^ -  -  0.008(2|av  +  yX-v)a3  + 


b(l  -  dv2) 


■b2032(i  -  nv2)2 


(107) 


Thrust  Unbalance 

Two-bladed  rotor.-  The  second -harmonic  variation  in  Gj.  for  a.  two- 
bladed  rotor  is 


AC, 


— —  =  Fourth  +  fifth  rows  of  equation  (53) 


(108) 


For  CDjj  =  toy  =  0  and  steady-state  conditions,  the  equation  for  the 
amplitude  may  he  simplified  to 

~  |_2a2I3s  ‘  (bl  +  |  w)M2s  +  \  Vv21  J  + 


(2b2l3s  +  alM2s)  *  j 


n  1/2 


(109) 
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Three -b laded  rotor The  third-harmonic  variation  in  Crp  for  a 
three -b laded  rotor  is  approximately 


2ACip 

- -  £  Sixth  +  seventh  rows  of  equation  (53)  (110) 

3a 


or  for  cj^  =  =  0  the  amplitude  is  approximately 


2ACrp 


V2lls 


(111) 


An  Independence-of -Blade-Element  Analysis  for  Hovering, 

Vertical  Ascent,  and  Convertaplane  Propeller  Condition 

The  use  of  the  relation  c^  =  a  sin  a  permits  a  considerable 
simplification  of  the  equations  resulting  from  the  assumption  of  the 
independence  of  blade  elements.  As  the  exact  propeller  solutions  of 
Betz,  Goldstein,  and  Theodorsen  are  not  applicable  to  a  lifting  rotor 
at  zero  or  small  advance  ratios,  a  simple  analysis  of  the  independence 
of  blade  elements  may  be  useful. 

From  momentum  considerations  the  thrust  dT  on  an  annulus  of  the 
rotor  disk  2itr  dr  is  related  to  the  induced  velocity  Vp  at  the  rotor 
element  by  the  expression 


dT 

4itpr  dr 


Vi(Vi  +  V  sin 


(112) 


But 


V-j^  +  V  sin  =  U  sin  0V 


(113) 


dT 

4npr  dr 


(u  sin  0v)(u  sin  0y 


-  V  sin 


(11*0 


Thus 


kb 
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The  thrust  of  the  annulus  is  also  equal  to  the  thrust  acting  on  the 
portions  of  the  blades  within  the  annulus  which  is 


dT  =  i  pbU2cc  j  cos  0y  dr 


(115) 


where 


c 2  ='  a  sin  =  a  (sin  0y  cos  0y  +  cos  0y  sin  0V)  (ll6) 


Thus 


dT  =  |  Pab (u  cos  0V)  [sin  0y(u  cos  0y)  +  cos  0y(u  sin  0y)]  c  dr  (117) 

Substituting  the  above  values  of  dT  in  equation  (ll4)  and  solving  for 
U  sin  0y 


U  sin  0y  /Vo  aba 


OR 


+  —r-  cos  ev 

16  t 


p 


va  abar  \2  aba~ 

+  cos  9vj  +  — g —  x  sin  0y 


(118) 


where 


va  = 


V  sin 
HR 


(119) 


jtR 


(120) 


Then  from  equation  (117) 


2C 

ab 


--It 


( U  sm  0y\ 

x  sin  0  +  I— -  cos  0lr 

v  V  HR  /  v 


arx  dx 


(121) 
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U  sin  0V 

where  the  value  of  -  at  x  is  given  by  equation  (ll8) 

HR 

Similarly,  from  blade-element  considerations 


1  /U  sin  0v\f~  /U  sin  0  \  "1 

V"--jR-  )[?  Sin  +  \  OR  7  c0s  6vJ^  dx  + 


sin  ar 


x  sin  0-ir  + 


U  Sin  ?v\  *~L  .2 


cos  a-rX^  dx 


(122) 


CLo  ao 

where  the  values  of  — -  are  obtained  from  a  plot  of  — : -  against 

sm  ar  sin  ar 

ar  for  the  blade  airfoil  at  values  of  ar  given  by  the  relation 


l[lfu  sin  0. 


0^  =  6^+  tan  x  — I - 


(123) 


If  it  is  necessary  to  take  into  account "the  rotation  of  the  slip¬ 
stream  for  large  rates  of  vertical  ascent  or  the  propeller  condition, 
this  may  be  accomplished  to  a  first  approximation  by  using  an  effec¬ 
tive  £2,  £2e,  in  every  case  where 


=  n  ^1  -  ^ 


(124) 


The  geometry  of  the  above  equations  is  exact  and  they  are  convenient 
for  graphical  or  numerical  integration  on  account  of  the  repetition  of 
factors. 

Neglecting  the  induced  radial  and  tangential  velocity  components, 
the  optimum  blade-angle  distribution  for  minimum  induced  power  and  a 
given  blade-chord  distribution  and  nondimens ional  axial  flight-path 

U  sin  0V 

velocity  va  may  be  obtained  by  setting  - ^ — -  equal  to  the  constant 

value  giving 
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sin  0V  = 


Xy  (Xy  -  Va)  X 

1  +  \ 

1  (  (\2  +  x2)[k2  -  (K  -  va)2]| 

k^2  +  x2) 

1 

w 

(\  ~  va)2*2  J 

(125) 


where 


and 


(126) 


(127) 


The  optimum  chord  dis trihut ion  for  a  given  desired  constant  value  of  c ^ 
along  the  blade  and  the  same  restrictions  is 


&*v(K  -  va) 

bcz^2  +  x2 


(128) 


i 

For  this  optimism  chord  distribution,  the  optimum  distribution  of  By 
reduces  to 


* 


sin  6V  = 


\ 


1  + 


a2Xv2  -  c +  x2) 


2.2 


C  7  x 


(129) 


.  For  calculations  where  the  flight-path  velocity  and  equilibrium 
value  of  Cfp  are  known  or  can  be  estimated,  the  following  procedure 
may  be  followed: 


(1)  Calculate  and  plot  the  radial  distribution  of 

(2)  'Calculate  the  effective  value  of 


CT~  "  CT 


=  V, 


f— V 
/ 

(<ft) 


Ciji  and  va  where 
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(3)  Calculate  the  approximate  value  of  A q  from  equation  (72) 
which  for  these  flight  conditions  reduces  to 

2CTe 

ab  “  03s  “  ^va2c 
(ct3c  "  ^v^s) 

(4)  Calculate  and  plot  the  radial  distribution  of  ev  =  Aq  + 
for  the  value  of  Aq  obtained  under  item  (3)  and  two  lower  values  at 
increments  of  several  degrees 

(5)  Calculate  and  plot  the  radial  distribution  of  .  U  sin  0v/fteR 

for  the  above  distribution  of  from  equation  (ll8)  using  Q  =  Qe 

throughout 

(6)  Calculate  and  plot  the  radial  distribution  of  the  integrand 
of  equation  (121)  for  the  three  values  of  Aq  and  graphically  or 

numerically  integrate  for  the  values  of  2CTe//ab  corresponding  to  the 
three  values  of  Aq. 

(7)  Obtain  the  correct  value  of  Aq  from  a  plot  of  CTg  against  Aq 

(8)  Calculate  and  plot  the  radial  distribution  of  the  integrand  of 
equation  (122)  for  the  three  values  of  Aq  and  graphically  or  numerically 

integrate  for  the  values  of'  SC^b  corresponding  to  the  three  values 

°f  Aq 

(9)  Obtain  the  equilibrium  value  of 
of  Aq  from  a  plot  of  CQg  against  Aq 

(10)  Calculate  the  equilibrium  value 


the  equilibrium  value 

°f  C«-C«e(lf)' 


sin  A, 


Comparison  of  Experimental  and  Calculated  Values  of  Parameters 

Table  6  shows  a  comparison  of  the  experimental  data  of  reference  2 
for  those  runs  where  CT  2  0.00545  with  the  values  calculated  by  the 

approximate  blade-element  equations  of  this  report.  The  blade-element 
lift-curve  slope  was  assumed  by  the  authors  to  have  been  a  =  6.5  from 
the  experimental  results  of  reference  6.  The  values  of  gq,  e±, 
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and  ^2  were  evaluated  for  the  points  =  0.0090,  0.0105,  and  0.0170 

at  c 1  =  0,  0*5,  and  1.0,  respectively,  from  figure  19  of  reference  6. 

The  exact  solutions  for  the  various  parameters  differ  from  the 
tabulated  approximate  solutions  by  a  negligible  amount  for  these  heli¬ 
copter  flight  conditions. 

A  consideration  of  the  results  presented  in  table  6  would 
indicate  that  much  of  the  remaining  discrepancy  between  experimental 
and  calculated  blade  angles  and  torque  coefficients  may  be  due  to  the 
neglect,  in  the  present  calculations,  of  the  effects  of  the  rotor  induced 
velocity  on  the  lift  and  drag  of  the  fuselage. 

It  may  be  noted  that  the  longitudinal  component  of  the  angle 
tan )  between  the  rotor  resultant  force  and  the  thrust  component 

\2ct/ 

normal  to  the  tip-path  plane  fh  very  small  for  all  these  helicopter 
flight  conditions  and  that  the  direction  of  the  resultant  is  inclined 
forward  for  those  flight  conditions  where  there  is  a  net  downflow  through 
the  rotor.  The  inclinations  of  the  tip-path  plane  to  the  horizontal  0X 
and  6y  are  also  small  angles  and,  consequently,  for  many  unaccelerated- 
flight  helicopter  calculations  the  rotor  resultant  force  can  be  assumed 
to  be  perpendicular  to  the  tip-path  plane  and  the  thrust  equal  to  the 
gross  weight  without  introducing  serious  errors. 


CONCLUDING  DISCUSSION 


Simple  relations  for  the  rotor  blade  angles  and  the  values  of  Cq, 
Cx,  and  Cy,  derived  upon  the  assumption  of  a  triangular  distribution 
of  blade-element  circulation  along  the  radius  and  a  sinusoidal  variation 
with  azimuth  angle  in  conjunction  with  a  linear  variation  of  profile 
drag  with  lift,  would  appear  to  be  useful  for  helicopter  and  converta- 
plane  performance  estimation  and  the  determination  of  the  equilibrium 
angle  of  attack  and  lateral  tilt  of  the  tip-path  plane. 

The  blade-element  equations,  based  upon  the  relation  that 
cl  -  a  sin  ar  =  a(sin  0V  cos  0V  +  cos  0v  sin  0V) ,  and  the  crnc  and 

ans  functions  of  the  blade-chord  and  blade-twist  distribution  afford 
a  reasonably  exact  and  concise  treatment  of  the  geometry  and  should  be 
useful  for  convertaplane  as  well  as  helicopter  calculations. 


s 
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The  use  of  the  empirical  relation  +  e\  sin  ar  +  ^2  cos  ar, 

o 

rather  than  the  usual  expression  that  c<i0  =  £>q  +  &]_ar  +  &2ar  >  consid¬ 
erably  simplifies  the  equations  for  the  in-plane  forces  and  moments  and 
presents  a  sufficiently  exact  solution  of  the  geometry  for  helicopter 
calculations.  *  , 

For  convertaplane  calculations,  the  approximation  that 
cd0  =  el  sin  +  €2  cos  “r  slices  an  exact  treatment  of  the  geometry 
and  should  be  a  sufficiently  accurate  expression  for  c^  at  the  larger 

advance  ratios  where  the  effects  of  the  profile  drag  become  of  less 
..relative  importance. 

The  larger  sources  of  the  remaining  errors  in  the  blade-element 
analysis  probably  have  the  following  order  of  importance  for  contemporary 
helicopters: 

(1)  The  neglect  of  the  effects  of  blade-element  stall  implied1*  in 

the  relation  that  c^  =  a  sin 

(2) N The  neglect  of  the  effects  of  blade  flexibility 

(3)  The  neglect  of  the  radial  variation  in  the  normal  component 

of  the  induced  velocity 

(4)  The  neglect  of  the  effects  of  compressibility  on  the  tip 

sections  of  the  advancing  blade. 


Georgia  Institute  of  Technology 
Atlanta,  Ga. ,  May  15,  1951 
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TABLE  1 

VALUES  OF  an  FOR  BLADES  WITH  LINEAR  TAPER 


^  o 

Interpolate  for  values  for  given  t;  an  =  — ; 

w  kR 

t  =  -  1;  c  =  cQ(l  +  t x)J 


t 

Q 

H 

Q 

O 

Q 

ro 

Q 

O 

°3/a0 

ak/a0 

X-,  =  0.15 

1  M 

0 

0.8500 

0.4888 

0.3322 

0.2499 

-1 

.3612 

.1566 

.0823 

.0499 

=  0.20 

0 

0.8000 

0.4800 

0.3307 

0.2496 

-1 

.3200 

.1493 

.0811 

.0497 
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TABLE'  2 

VALUES  OF  crnc  FOR  BLADES  WITH  LINEAR  TAPER,  LINEAR 
TWIST,  AND  x±  =  0.15 


I 


Interpolate  for  values  for  given  t  first  and  then  for  values  for 


given  0-^;  reference  station  for  Aq  at  x  =  0; 


_  ctjp 


l; 


c  =  cQ(l  +  tx); 


9, 


=  0-.X 


a  -  - o. 
u  JtR 


el 

(deg) 

CTlc/a0 

ff2c/a0 

CT3c/ao 

a4c/ 

ao 

t  =  0 

m 

t  =  0 

m 

ns 

t  =  0 

t  =  -1 

0 

0.8500 

0.3612 

0.4888 

0.1566 

0.3322 

0.0823 

0.2499 

0.0499 

-4 

.84-92 

.3611 

.4882 

.1565 

.3317 

.0822 

.2495 

.0498 

-8 

.8468 

.3604 

.4864 

.1561 

.3303 

.0820 

.2483 

.0497 

-12 

.8427 

.3594 

.4833 

.1555 

.3278 

.0816 

.2462 

.0494 

-16 

.8371 

.3580 

.4791 

.1546 

.3244 

.0810 

.2434 

.0490 

-20 

.8299 

.3562 

.4737 

.1536 

.3201 

.0803 

.2398 

.0485 

-24 

.8211 

.3541 

.4671 

.1522 

.3148 

.0794 

.2354 

.0478 

-28 

.8108 

.3515 

.4594 

.1507 

.3087 

.0784 

.2303 

.0471 
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TABLE  3 

VALUES  OF  cxns  FOR  BLADES  WITH  LINEAR  TAPER,  LINEAR 
TWIST,  AND  x1  =  0.15 

Interpolate  for  values  for  given  t  first  and  then  for  values  for 

0 

given  0-, ;  reference  station  for  An  at  x  =  0;  crn  =  — £■; 

x  u  v  jtR 

t  =  -  1;  c  =  c0(l  +  tx).;.  6t  =  S-lxJ 


el 

CTls/a0  ' 

<J2s/cr0 

a3s/a0 

°4s/a0 

(deg) 

t  =  0 

t  =  -1 

t  =  0 

t  =  -1 

t  =  0 

t  =■  -1 

t  =  0 

t  =  -1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-4 

-.0341 

-.0109 

-.0232 

-.0057 

-.0174 

-.0035 

-.0139 

-.0023 

-8 

-.0681 

-.0219 

-.0463 

-.0115 

-.0348 

-.0070 

-.0279 

-.0047 

-12 

-.1020 

-.0327 

-.0693 

-.0172 

-.0521 

-.0104 

-.0417 

-.0070 

-16 

-.1356 

-.0435 

-.0920 

-.0229 

-.0692 

-.0138 

-.0553 

-.0092 

-20 

-.1699 

-.0531 

-.1145 

-.0284 

-.0860 

-.0173 

- . 0688 

-.0115 

-24 

-.2017 

-.0650 

-.1367 

-.0341 

-.1026 

-.0206 

-.0820 

-.0137 

-28 

-.23^0 

-.0756 

-.1585 

-.0396 

-.1189 

-.0239 

-.0950 

-.0159 
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TABLE  4 


VALUES  OF  ffnc  FOR  BLADES  WITH  LINEAR  TAPER,  HELICAL 


TWIST,  AND  =  0.20 


Interpolate  for  values  for  given  t  first  and  then  for  values  of 

r* 


given  reference  station  for  AQ  at  blade  tip;  cr 


•o. 


t  = 


'tip 


-  1}  c  =  c  (1  +  tx);  0- 


tan 


_1(x  tan  Qt} 


0  JtR'* 


0T 

alc/a0 

Q 

ro 

0 

Q 

O 

o3c/°0 

a4c/°o 

(deg) 

t  =  0 

B 

t  =  0 

t  =  -1 

B 

t  =  0 

t  =  -1 

0 

0.8000 

0.3200 

0.4800 

0.1493 

0.3307 

0.0810 

0.2496 

0.0497 

-4 

.7906 

.3144 

.4762 

.1474 

.3287 

.0803 

.2484 

.0493 

-8 

.7654 

.3002 

.4651 

.1419 

•  3233 

.0782 

.2451 

.0483 

-12 

.7305 

.28o4 

.4500 

.1351 

•  3149 

.0751 

.2398 

.0468 

-16 

.6907 

.2594 

.4313 

.1270 

.3042 

.0714 

.2328 

.0448 

-20 

.6489 

.2385 

.4io4 

.1185 

.2919 

.0674 

.2244 

.0427 

-24 

.6065 

.2184 

.3882 

.1100 

.2782 

.0632 

.2151 

.o4o4 

-28 

•5645 

.1994 

.3651 

.1016 

.2635 

.0590 

.2046 

.0379 

-32 

.5231 

.1815 

.3416 

.0934 

.2481 

.0547 

.1935 

.0354 
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TABLE  5 

VALUES  OF  ong  FOR  BLADES  WITH  LINEAR  TAPER,  HELICAL 
TWIST,  AND  xx  =  0.20 


1“ 


^Interpolate  for  values  for  given  t 
given  0rp;  reference  station  for 


first  and  then  for  values  for 
A0  at  blade  tip;  aQ  = 


t 

ctip 

C° 

H 

V.  • 

0 

II 

c0(l  + 

tx) ;  @t 

=  tan  ^ 

(x  tan  6t)] 

0rp 

cls/a0 

02s/aO 

°3s/c0 

CT4s/ao 

(deg) 

t  =  0 

t  =  -1 

t  =  0 

t  =  -1 

t  =  0 

t  =  -1 

t  =  0 

t  =  -1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-4 

-.1106 

-.0553 

-0553 

-.0220 

-.0333 

-.0103 

-.0230 

-.0056 

-8 

-.2121 

-.1045 

-.1076 

-.0422 

-.0654 

-.0200 

-.0454 

-.0110 

-12 

-.3005 

-.1452 

-.1553 

-.0596 

-.0957 

-.0287 

- . 0669 

-.0160 

-16 

-.3761 

-.1780 

-.1981 

-.0744 

-.1237 

-.0364 

-.0872 

-.0205 

-20 

- . 44o5 

-.2044 

-.2365 

-.0872 

-.1494 

-.0431 

-.1062 

-.0245 

-24 

-.4956 

-.2256 

-.2701 

-.0972 

-.1728 

-.0490 

-.1239 

-.0281 

-28 

-.5430 

-.2428 

- . 3002 

-.1060 

-.1941 

-.0540 

-.l4oi 

-.0313 

-32 

-.5838 

-.2570 

-.3269 

-.1134 

-.2134 

-.0584 

-.1551 

-.0342 
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COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  VALUES  OF  PARAMETERS  FOR 


c3  sq  ft  added  to  fuselage  drag  area  f  from  reference  7  to 
allow  for  rotor-hub,  blade-shank,  counter- torque -rotor,  and  engine-cooling  drag. 

^Includes  correction  for  tip  stall. 

Mechanical  input  subtracted. 


COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  VALUES  OF  PARAMETERS  FOR 


ti 


H  bo 
ri  a 
•H 

O 

V 


u 

o 


Mechanical  input  subtracted. 
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TABLE  7  -  Concluded 


VALUES  OF  FOR  GIVEN  VALUES  OF  A*  =  M^'2  “  3Mv 


UR 


T 


V  sin  cl,.  2  -  3m/ 

AND  =  - - — —ll - — - Concluded 


fiR 


'T 
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Z(T,-Q,fl,i|r) 

i 


Figure  1.-  Tip-path  plane  or  axes  of  virtual  rotation. 
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calculations  and  determination  of  equilibrium  angle  of  calculations  and  determination  of  equilibrium  angle  of 

attack  and  lateral  tilt  of  the  tip -path  plane.  More  attack  and  lateral  tilt  of  the  tip -path  plane.  More 

exact  blade -element  equations  are  then  derived  for  exact  blade -element  equations  are  then  derived  for 

the  relations  between  thrust  and  flight -path  velocity  the  relations  between  thrust  and  flight -path  velocity 

components  and  equilibrium  blade  angles,  torque,  components  and  equilibrium  blade  angles,  torque, 
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